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The s = | Ising spin chain with uniform nearest-neighbor coupling, quadratic single-site potential, 
and magnetic field is shown to be equivalent to a system of 17 species of particles with internal 
structure. The same set of particles (with different energies) is shown to generate the spectrum of 
the s = | Ising chain with dimerized nearest-neighbor coupling. The particles are free of interaction 
energies even at high densities. The mutual exclusion statistics of particles from all species is 
determined by their internal structure and encoded in a generalized Pauli principle. The exact 
statistical mechanical analysis can be performed for thermodynamically open or closed systems and 
with arbitrary energies assigned to all particle species. Special circumstances make it possible to 
merge two or more species into a single species. All traits that distinguish the original species 
become ignorable. The particles from the merged species are effectively indistinguishable and obey 
modified exclusion statistics. Different mergers may yield the same endproduct, implying that 
the inverse process (splitting any species into subspecies) is not unique. In a macroscopic system 
of two merged species at thermal equilibrium, the concentrations of the original species satisfy a 
functional relation governed by their mutual statistical interaction. That relation is derivable from 
an extremum principle. In the Ising context the system is open and the particle energies depend on 
the Hamiltonian parameters. Simple models of polymerization and solitonic paramagnetism each 
represent a closed system of two species that can transform into each other. Here they represent 
distinguishable traits with different energies of the same physical particle. 



I. INTRODUCTION 

The idea that matter consists of particles, though an- 
cient, has remained controversial until the dawn of the 
twentieth century. All known particles are classified into 
species hierarchically, e.g. at the level of quarks, at the 
level of baryons, mesons, and leptons, at the level of 
atoms, and at the level of molecules. The indistinguisha- 
bility of particles from the same species - a key concept 
of statistical mechanics that naturally arose from symme- 
try requirements in quantum mechanics - has clear-cut 
consequences in the microscopic and macroscopic realms 
of physical reality [TJ [5] . The universality of this concept 
is supported by overwhelming evidence. 

The notion of particle species in the context of indistin- 
guishability is subtle. Two electrons may be distinguish- 
able by their spin orientation, their momenta, or their 
angular momenta. Different isotopes distinguish atoms 
of the same element, and different isomers distinguish 
molecules of the same composition. All distinguishable 
traits are observable attributes. Particles that share all 
traits are identical. However, particles need not be identi- 
cal to belong to the same species. Not all distinguishable 
traits are relevant in all contexts. 

For the purpose of this study we use motifs to charac- 
terize distinguishable traits of particles. Particle species 
with a single motif are classified at the most fundamental 
level. At levels of classification that ignore certain traits, 
particles belonging to one species are characterized by 
one of several motifs. Inevitably, the particle content 
(numbers of particles from each species) of many-body 
states is of a different mix at different levels of classifi- 
cation. However, the statistical mechanical analysis of a 
system of particles must produce consistent results at all 



levels of classification. Understanding these consistency 
criteria is one of our goals. 

Aiming to hide all distinguishable traits that do not 
affect the physical quantities under investigation mini- 
mizes the number of species and assigns multiple motifs 
to some or all of them. It is one optimizing factor for 
the analysis. The second optimizing factor selects par- 
ticle species such as to weaken or eliminate interaction 
energies between particles. The two operate in the same 
arena but have different ends and are subject to different 
constraints. 

In two previous papers we established a scheme of clas- 
sifying statistically interacting particles with shapes into 
species according to structure and into categories accord- 
ing to (primitive) functions. We investigated the crite- 
ria that make the particles free of interaction energies 
[31 13] . In this work we demonstrate how multiple species 
of statistically interacting particles are merged or split 
according to need and economy in a way that produces 
consistent and exact results at all levels of classification. 
We again select a particular many-body system to de- 
velop an idea of much wider scope. 

We consider the spin-| Ising chain with nearest- 
neighbor coupling, a quadratic single-site potential, a 
magnetic field, and periodic boundary conditions: 
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where Sf = ±| , ±™. It is sufficiently complex to fa- 
cilitate multiple levels of classification for the particles 
assembled by the spin coupling. For the product eigen- 
states we use the notation \<Jio~2 ■ ■ ■ <jn)p, where the pe- 
riodicity p denotes the number of distinct vectors related 
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via translations. To the site variables oi we assign sym- 
bols as follows: (+§, +|, -|, -§) = (ff, t, t, |1). 

There exist four zero-temperature phases at h = in 
sectors of the (J, Z?)-plane: 

: | ftt •■•>!. I 444-")i at J<0,£>> - J; 
$f : at J<0,£><-J; 

$a : I tit •••>2 at J> 0,D > J; 

|flf->2 at J>0,£>< J. (2) 

At /i ^ the number of phases increases to ten [5]. 
All of them are realized for J > 0. Their ranges in the 
reduced parameter plane (D/ J, h/J) are shown in Fig.JT] 
The magnetization increases from zero to saturation in 
one step if D/J < and in three steps if D/J > 0, The 
plateau phases are at one third and two thirds of the sat- 
uration magnetization. Four phases are non-degenerate 
and six are twofold degenerate. When h increases from 
zero at Dj J > 1 the translational symmetry increases at 
the first transition, decreases at the second, and increases 
again at the third. In the limit D — > oo the second term 
in ([I]) reduces the Hilbert space of finite-energy states to 
that of the s — | Ising spin chain. 
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FIG. 1: T = phases of model system {I}. 

In Sec. [n] we introduce, for the s = | Ising chain, one 
set of 17 species of particles that remain free of interac- 
tion energies even at high densities. Two species do the 
same trick in the limit D — > oo. In Sec. Mil we formu- 
late criteria for merging two or more species into one via 



combinatorial analysis. In Sec. IV we discuss mergers of 
species in the context of the statistical mechanical anal- 
ysis with applications to open systems. Closed systems 
of particles with distinguishable traits are discussed in 



Sec. [Vj In Sec. VI the major themes are recapitulated 
and their threads tied together. 



II. FLOATING PARTICLES 

All A N product eigenstates of , arrays of er; , are now 
interpreted as strings of particles and elements of pseudo- 
vacuum, represented by motifs that comprise between 
two and four consecutive site variables in distinctive pat- 
terns. Motifs interlink along the chain in one shared site 
variable. 

As in previous applications to Ising spin chains with 
s = | , 1 [21 IH E] , the coupling of the oi causes the as- 
sembly of particles from a set of species that are free of 
interaction energies and generate the entire spectrum of 
@. The pseudo -vacuum of choice in this work will be 
the twofold Neel state | f|-|!f|- ■ ■ • )i- It coincides with the 
physical vacuum in left center of the parameter space 
depicted in Fig. [TJ 

The search for particles that remain free even at high 
density is facilitated by the taxonomy introduced in 
Ref. [5] using familiar and well-tested notions of frac- 
tional exclusion statistics [3[S]. For the multiplicity func- 
tions W({N m }), which expresses the number of product 
states that contain specific numbers {N m } of particles 
from all species, we use the same template with no need 
for adjustments: 
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where n pv is the multiplicity of the pseudo-vacuum, the 
A m are capacity constants, the a m are size constants, and 
the g m m' are statistical interaction coefficients. Relation 
( 3b ) embodies the generalized Pauli principle as proposed 
by Haldane [7 . The integer d m counts the number of 
open slots for particles of species m in the presence of 
N m i particles from any species ml '. 

The absence of any binding energies, which allows us 
to view the particles as floating objects, indistinguishable 
except for species, is reflected in the energy expression, 



E({N m }) = E, 



pv 



M 
m— 1 



(4) 



for product states with particle content {N m } and parti- 
cle energies e m measured from the pseudo-vacuum. The 
specifications to be used in ([3| and Q of the M = 17 
species assembled by the Ising spin coupling in are 
compiled in Table [T] 
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TABLE I: Specifications of M = 17 particles that generate the spectrum of |T]) from the state (n pv = 2), | 444 ■ ■ ■ )2- motif, 
species, energy (relative to vacuum at h = 0), spin, capacity constant, size constant, and statistical interaction coefficients. 
Segments of I vacuum bonds 44,44 have energy E = ^£(8D — 9J). For h 7^ the entries of e m must be amended by —s m h. 
Species m = 1, 2 are compacts, m = 3, . . . , 14 hosts, and m = 15, 16, 17 hybrids. 
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To determine which categories of particles are repre- 
sented in this set we recall [3] that compacts and hosts 
float in segments of pseudo-vacuum, whereas tags are 
located inside hosts, and hybrids are tags with hosting 
capability or tags that facilitate coexistence with other 
tags inside the same host. Hence the first two species 
are compacts, and the last three are hybrids. All other 
species are hosts. 

The hybrids alone generate 2 N product states, each one 
jammed with particles. The remaining 4 N — 2 N states - 
the vast majority in a macroscopic system - contain at 
least one compact or host or element of pseudo- vacuum. 
Four of the hosts are characterized by two motifs each, 
an attribute encountered previously [21 HI [S] , which does 
not cause a counting ambiguity because at most one of 
the two motifs can be placed into any particular slot. 

We also consider the spectrum of finite-energy states in 
the 2 w -dimensional Hilbert space pertaining to the limit 
D oo. Two sets of M = 2 species generate that spec- 
trum from different pseudo- vacua. The combinatorics of 
both sets are governed by multiplicity expression ([3| with 
specifications as compiled in Tables [TT] and |III| One set 
features two compacts, the other a host and a tag. Both 
sets were previously encountered as a special case in a 
different context [J]. 



TABLE II: Specifications of M = 2 species of particles excited 
from from the Neel state (n pv — 2) | tltl ' ' ' )i- Segments of 
I vacuum elements, tl,lt, have energy —U/4. At h 7^ the 
entries of e m must be amended by —s m h. 



motif 


cat. 


m 






A m 


a m 


Qmm 1 


+ 




tt 


comp. 


+ 


J 

2 


+ 1 


iV-l 

2 


1 


+ 


1 
2 


1 
2 




comp. 




J 
2 


1 
2 


JV-l 
2 


1 




1 
2 


1 

2 



TABLE III: Specifications of M — 2 species of particles ex- 
cited from the spin-polarized state for (n pv = 1) | tt ' ' ' )i- 
Segments of £ vacuum elements, ft, have energy U/A. At 
h 7^ the entries of e m must be amended by —s m h. 
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Particles with the same motifs as in Table UJ but dif- 
ferent energies are relevant for the s = \ Ising model 
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with dimerized nearest-neighbor coupling as described in 
Appendix [A] 



III. SPECIES MERGED OR SPLIT 

In the following we investigate the circumstances un- 
der which particle species with multiple motifs can be 
produced via mergers from particles with a single motif. 
It is also possible, albeit not without ambiguity, to split 
a species with multiple motifs into multiple one-motif 
species 0[TO]- The key criterion is that, at each level of 
description, particles from all species can be treated as 
indistinguishable, either because they are identical (sin- 
gle motif) or because their distinguishable traits (mul- 
tiple motifs) are not relevant in the given context 
We shall see that the criteria are more stringent on the 
level of combinatorial analysis (for systems of all sizes) 
than on the level of statistical mechanical analysis (for 
macroscopic systems). 



are identical if condition (5e) is satisfied. We are thus 
left to prove the equation 



d + N - 1 
N 




with 



x d 2 + N - Nl -l\ 
\ No-Ni ) 



d + N - 1 = D + g 00 - 1 - (.9oo - 1)JV , (9) 



d 1 + Ni-l = D 1 +g n -l- g 12 N 
+ (.9i2-5n + l)^i, 



(10) 



d 2 + No - JVi - 1 = D 2 + g 22 - 1 - 321^0 

+ (fl2i-flaa + l)(JVo-JVi) J (H) 



A. Rules 

Consider a set of M species with all specifications that 
go into the multiplicity expression ^ and the energy 
expression Q given. If, say, the first two species m = 1, 2 
are to be merged into a new species m — 0, this is possible 
if their specifications satisfy Anghel's rules [5] and are 
assigned values as follows: 



(i = £2 = eoj 
a i = a 2 = ao, 
Ai +A 2 = A , 

.911 + 521 = .912 + 522 = 5oo, 

5ml = 5m2 = 5mo, m = 3, . . . , M, 



5lr, 



52r. 



50r, 



m = 3, 



,M. 



(5a) 
(5b) 
(5c) 
(5d) 
(5e) 
(5f) 



Anghel derived these rules in an effort to ensure consis- 
tency in the statistical mechanical analysis of a system 



at different levels of coarse-graining [10] . Rules ( 5c ) and 



(5d) have to be made more stringent for the combina- 



torial analysis of systems of any size. The multiplicity 
expressions (3a) and the energy expressions Q of the 
old and new sets of species must satisfy the relations 

W(N ,N 3 ,...) = W{Ni,N 2 ,N 3) ...), (6) 

{Ni,N 2 } 

E(N ,N 3 ,...)= E(Ni,N 2 , N 3 , . . .), (7) 

{ATi,iV2} 

respectively, where the sum is constrained to pairs with 
Ni+N 2 = No as implied by the definition of the merging 
operation. Relation (JtJ) only depends on rule (5a). 

In proving relation ^ we consider all factors of (3a) 
separately. The prefactors left and right are identical if 
condition (5b) holds. All binomial factors with m > 3 



D, 



^2g mm >N m ,, m = 0,l,2. (12) 

m'— 3 



We note that D\ + D 2 — Do if conditions ( 5c ) and fl5 



are satisfied. Next we enforce condition ( 5d I by means 
of the parametrization 



1 



1 



5ii = 2-9oo + u, 512 = 2-9oo + «, 

521 = ^5oo - u, 522 = ^.9oo - v. (13) 

We have identified two scenarios in which {[sj) is proven 
to hold. Both necessitate a tightening of condition (5d) 



and one requires that condition ( 5c ) be tightened as well 



Type 1 mergers require that u = v+1, thus eliminating 
the last term in both (10) and (11) and reducing Q to 
the identity 




m + n 
k 



(14) 



Type 2 mergers require that u — v — and A\ = A 2 , 
implying D\ = D 2 and reducing ^ to the identity 




2m + k + 1 



(15) 



k — i I \ k 
The amendments to rules |5]) in generic form thus read 

3il = 512 ±1, 522=521 + 1 (type 1), (16a) 
5n = 5i2 = .921 = .922, Ax = A 2 (type 2). (16b) 

Anghel's rules ^ are determinate when species are 
merged but some are indeterminate - even with amend- 
ments ( 16 ) - without implied assumptions or additional 



input when species are split [9TU1|. 
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B. Applications 

We first apply Anghel's rules ([5| with amendments 
(161 to the pairs of species listed in Tables [TT III The 
two compacts can be merged if e + = e_, which is the 
case at h = and arbitrary values of J. The same is true 
for the host and the tag if e# = ex, which is the case at 
J = and arbitrary values of h. 

Both mergers, the first being of type 2 and the sec- 
ond of type 1, produce one species of lattice fermions 
with specifications Aq = N — 1, goo = 1, (Xq = 1, sig- 
nalling fermionic statistics. The distinguishable traits of 
the original species are erased from the specifications A m , 
gmm' that go into the statistical mechanical analysis. 

In both cases it is possible to replace the two motifs of 
the original species by a single motif representing float- 
ing lattice fermions. In the merger of host and tag, a 
(fermionic) particle (•) is assigned to any down-spin and 
a hole (o) to any up-spin. The mapping is one-on-one 
with a non-degenerate pseudo-vacuum. In the merger of 
the two compacts a particle is assigned to any aligned 
bond and a hole to any anti-aligned bond. The mapping 
is two-on-one with a twofold pseudo- vacuum |12j . 

If our goal is to investigate the statistical mechanics of 
the 17 species of particles listed in Table [I] (via methods 
explained in Sec. IV) we can simplify the analysis by first 



merging all species whose distinguishable traits are not 
relevant in a given context. We note that five species 
already have two motifs each, which means that some 
of their distinguishable traits happen to be redundant in 
the context of model system [T3] . 

For the special case h = we find that all twelve species 
with one motif can be sorted into pairs to be merged se- 
quentially on the basis of Anghel's rules §5§ and amend- 
ment ( 16b ). Not all sequences are permissible but all per- 



missible sequences yield the same results. One allowed 
sequence of mergers is the following: 

15 k 16 -> 15, 1 & 2 -> I, 5 & 6 -> 5, 

3 & 4 -> 3, 9 & 10 -> 9, 11 & 12 ^ fl. (17) 



Species 5 and 6, for example, do not satisfy (5e) be- 
fore the two preceding mergers have been carried out. 
The six mergers produce 11 species with two motifs each 
and specifications as compiled in Tabic |IV| None of the 
mergers can be undone on the basis of the specifications 
£m,A m ,a m ,g mm < from Table IV and the rules ([5]), ( |16b ) 
alone. At least some of the contents of Table U will have 
to be worked out from scratch via combinatorial analysis 
of the given motifs. 



IV. OPEN SYSTEMS 

In the statistical mechanical analysis of the particles 
listed in Tables [T| |IV| the grandcanonical ensemble is the 
natural choice. The average populations of all species are 
controllable via chemical potential. Macroscopic systems 



TABLE IV: Specifications of M = 11 particles produced via 
six mergers from those listed in Table [I] At h = the spin s m 
is an ignorable trait. 
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that are closed for individual species or combination of 
species can then be constructed by a switch of indepen- 
dent variables on the level of thermodynamic functions 
(Sec. [v]). In the Ising context the system is open and the 
chemical potentials are absorbed in the particle energies. 



A. Method 



Our starting point is the expression, 
Z=J2 ^({^ m })exp(-]T 

{N m } \ m 



(18) 



for the grand partition function of a system of statis- 
tically interacting particles with given energies e m and 
given specifications A m , a m , g mm ', n pv that go into the 
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multiplicity function Wu's analysis of (18) for a 

generic situation [14j produced the general result 



1 + w Tl 



where the (real, positive) w m are solutions of 

' 1 + uw 



= ln(l + w m ) - 22 g m > m In 



(19) 



(20) 



The average numbers of particles from each species are 
inferred from 



w m {N m ) + } j gmm>{N m <) = A r , 



(21) 



The entropy derived from ( 18 ) can be expressed as a func- 
tion of the (N m ) alone: 

m 

- (iV m > ln(iV TO ) - (Y m ) hx(Y m )] , (22a) 

^rn A m ^ ^ Qmra' ^ra' ■ (22b) 
ni ' 

Merging species in the statistical mechanical analysis 
by erasing distinguishable traits of no relevance simplifies 
the calculation without sacrificing rigor. If the specifica- 
tions of any two species, say m = 1,2, obey Anghcl's 
rules (|5|, the consequences are that 



Wl = w 2 = w a , 
(JVi) + (N 2 ) = (N ) 



(23a) 
(23b) 



satisfy Eqs. ( 19 )-( 21 ) with 



and that species m — 0, 3, . 
the modified A m , g mm > ■ 

Interestingly, the additional restrictions imposed by 
( 16 1 are ignorable. Implementing the merger in the sta- 



tistical mechanical framework is less restrictive than in 
the combinatorial framework. The former is good for 
macroscopic systems, the latter for systems of all sizes. 
Type 1 and type 2 mergers leave distinctive signatures 
even in the statistical mechanical analysis, nevertheless. 

The entropy of mixing of two species of particles with 
shapes does not necessarily have its maximum when both 
species are present in equal numbers or equal concentra- 
tions. This fact is evident when we compare the entropy 
expressions ( [22] ) before and after a merger. Both can 
be calculated independently from the A m ,g mm i before 
and after any merger carried out in the framework of the 
combinatorial analysis (Sec. |TT| ) . Consistency of the two 
entropy expressions implies a functional relation between 
(Ni), (N2) in addition to (23b). It can be derived from 
the extremum principle, 



d 



S((N 1 ) ) (N )-(N 1 ) ) (N 3 ),...) 



= 0. (24) 



B. Two species 

In Fig. [2]_we show contour plots of the entropy land- 
scapes S((N+), (AL)) and S((N H ),(N T )) pertaining to 



the two pairs of particle species from Tables [TT] and III 
We use the notation S = S/N for the entropy per site and 
(N m ) = (N m ) /N for the average population densities. In 
generic open systems the (N m ) are functions of T, e + , e_ 
or T,€H,eT- In the Ising context, more specifically, the 
particle energies are functions of the Hamiltonian param- 
eters J, h. 

Both landscapes have triangular shape [TS]. Panel 
(a) tells us that each soliton species is able to produce 
entropy on its own (through mixing with elements of 
pseudo-vacuum) and yet higher entropy through mixing 
with each other at moderate densities. However, when 
the system is crowded with solitons, the entropy is low 
no matter what the ratio of the two species is. Spin- 
up and spin-down solitons are segregated at the highest 
density. 

From panel (b) we learn that only the hosts can pro- 
duce entropy on their own, not the tags, which exist in- 
side hosts in uniform arrays. Hosts with tags with tags 
inside come in many different sizes and thus produce yet 
higher entropy. The entropy stays high even when the 
system is close-packed with particles from both species. 




0.6 0.8 



(No) 



FIG. 2: Entropy per site, S/ks, versus average popula tion 
densities of the two species from (a) Tableland (b) Table|ITTl 
The contours are at (a) 0.070^, (b) 0.069^, £ = !,..., 9. The 
dashed lines represent relations (1261). 
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Implementing the extremum principle ( 24 ) yields iden- 
tical entropy expressions, 

S/k B = -(N Q ) ln(N ) - (l - (No)) In (l - (N )) , (25) 

reflecting fermionic exclusion statistics, but different 
functional relations between the average population den- 
sities of the merged species: 



(Nh) 



(N-), 



(N T ) - (N T ). 



(26a) 
(26b) 



These relations are independent of the particle ener- 
gies. They are solely governed by the shapes of the 
particles and the way they interlink. Adding fermions 
with hidden traits of spin-up and spin down solitons pro- 
duces a fifty-fifty mix at all densities. In strong contrast, 
fermions with hidden host-tag traits are mostly hosts at 
low density and mostly tags at high density. 



C. Seventeen species 



The statistical mechanical analysis of ([!]) at h ^ 



cannot take advantage of any mergers. Coupled Eqs. ( 20 ) 
and (21) must be solved for m = 1, . . . , 17. Expression 



( 19 1 for the grand partition function can, nevertheless, 



be written as a function of a single w„ 



[(l + Wl )e- 



18Kj+6H 



N 



(27) 



where w\ along with W2, ■ ■ ■ , 1017 are functions of Kj = 
J/4k B T, K D = D/4k B T, H = h/2k B T. 

In Fig. [3] we show the dependence of various expecta- 
tion values on the magnetic field along two parallel paths 
in the parameter space of Fig. IT] The aver age population 
densities (N m ) inferred from ( |20| and (21) determine the 
magnetization via 



My 



,(N m ) 



(28) 



and the entropy via ( 22 ) 



In each case there are two magnetization plateaux, 
softened by thermal fluctuations (shown for one case). 
The entropy shoots up at all three phase boundaries. 
Only the third phase boundary is common to both cases. 
In each phase region except one the population density 
of one species is high. Here the ground state is a solid of 
that species of particles. In one phase region, the ground 
state is the pseudo-vacuum, which, by definition, is de- 
void of particles. 

The ground-state degeneracy is huge at all phase 
boundaries crossed by the two paths. Particles or ele- 
ments of pseudo- vacuum with high densities in adjacent 
phase regions behave like miscible liquids with each spike 
in the dashed lines representing an entropy of mixing as 
discussed in Ref. 3 for a different model. 





1 2 3 4 5 6 
h/J 



FIG. 3: Entropy S/ks per site and average population den- 
sities (N m ) for selected species versus h/J at ksT/J = 0.1 
and (a) D/J = 0.6, (b) D/J = 1.4. The inset shows the 
magnetization M z per site versus h/J for case (a). 



Many more ground states than the ten phases shown 
in Fig. [l] can be stabilized by assigning different energies 
e m to the particles in Table [IJ This includes ground states 
that represent solids of more than one species in ordered 
or disordered configurations. 

If we assign the lowest energy (per bond) to particle 
m = 7, for example we produce a ferrimagnetic state 
with periodicity p = 4: | "frt-IJ-t 1 • • )4- If we assign the 
same lowest energy per bond also to particle m = 8 we 
produce a ground state that is antifcrrromagnetic but has 
a Langevin paramagnet embedded in it: | -ft a JJ. a ■ ■ ■ ) 
with cr =y, j in a random sequence. For each scenario 
thus designed the statistical mechanical analysis can be 
performed exactly. 



D. Eleven species 



At h — we can take advantage of two layers of simpli- 



we have merged six 
, thus reducing that 



fications in the analysis. In Sec. KIT 
pairs of the M — 17 species (Tablell 
number to M = 11 (Table IV), each with two motifs 
The simplifications do not stop here 
Inspection of the coefficients g m m 



in Table IV shows 



that five pairs of species almost qualify for further merg- 
ers. The pairs 7 & 8, 9 & 11 only violate condition |5dl, 
while pairs 3 & 5, 11 & 13, 13 & 14 violate also condi- 
tion (5a I. Even though these five pairs of species with 

9mm 9mm' 9m'm f — 9m'm ^) ^m ^-m' ^mm 1 

cannot be merged, the intact conditions from ^ lead 
to further simplifications in Eqs. (20) and plj ) for those 
pairs: 



1 + w r . 



{N„ 



exp 



A., 



(29) 



With all these simplifications in place, Eqs. (20) in 
exponential form can be rendered more compactly: 



l + wi 

1 + Wf 5 



— e 8Kj+8K D 
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W 15 _ „4.Kj 



1 + wn 



w sWf 5 = el2Kj wiw$ = el2Kj 

wir 1 + w$ 

1 + w 8 



,12if j 



1 + W3 



w s w 14 iK w 13 12K 

, =6 — C 

w 7 (l + w 7 )wn l + wu 



1. 



1 + tm _ ex.; 1 + _ 

W 13 Wg 

(1 + w 7 )(l + wg) = gl2ifj 



(30) 



and the grand partition function again reduces to (27), 
now for wi = u>j and H = 0. 

For the discussion of our results we consider a circular 
path of unit radius in the (J, D) plane. In Fig. [4] we 
show the variation of the entropy ( 22 ) along that path at 



various various temperatures, the same arbitrary units 
are used for J, D, and kgT. Naturally, this quantity is 
sensitive any border crossings of the path between the 
T = phases However, that sensitivity has some 
unusual features. 




FIG. 4: Entropy S/ks per site along a circular path, 
J = cos8, D = sin (9 in the (J,D) plane at various values 
of ksT in arbitrary energy units. The boundaries between 
the T — phases |2| are marked by vertical bars. 
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0.0 0.2 0.4 0.6 0.8 

6/7T 



1.0 



FIG. 5: Entropy S/ks per site and population densities (JV m ) 
along a circular path, J = cos#, D = sin#, (a) at fesT = 0.5 
(full circle) and (b) at ksT — 0.125 (half circle) in arbitrary 
energy units. 



The ground-state degeneracy is large at the $ a — $/ 
and $^4 — $i? boundaries but small at the <E> Q — 5>a and 
$/ — boundaries. At fcsT = 0.5 the entropy S has 
smooth maxima of moderate width at three borders and 
a broad minimum at the fourth. Surprisingy, an entropy 
minimum is located at the $ a — $/ border associated 
with a highly degenerate ground state. 

Upon lowering the temperature, these extrema evolve 
in different ways. The entropy maximum at the <&a — &f 
border remains high and sharpens. The entropy mini- 
mum at the $ a — $/ border turns into a maximum of 
equal height but of larger width. The maxima at the 
$>A — &a and $/ — $f borders shrink quickly. 

Only three of the eleven species, the compact 1 and the 
two hybrids 15 and 17, are present at high densities in 
the (J, D) plane. Their variation along the circular path 
is shown in Fig. [5] in combination with entropy data. 

The behavior near the $f — border is of the kind 
already discussed in the context of Fig. [3] compacts T and 
elements of pseudo- vacuum act as miscible liquids. The 
contrasting behavior of the data near the &a — ®a and 
$/ — $f borders can be explained by the observations 
(discussed in Ref. [3] for a different model) that hybrids 
17 and elements of pseudo-vacuum or hybrids 15 and 
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compacts 1 act as immiscible liquids. The population 
densities vary abruptly and the entropy of mixing is very 
small at low T. 

What is unusual about the behavior near the <E> Q — $/ 
border is that the energies of hybrids 15 and hybrids 17 
happen to cross at a very shallow angle. Both species 
maintain high population densities across both phases 
except at very low T. The mixing is not significantly 
enhanced near the border at moderate T. Hence the 
absence of a spiked entropy. 



V. CLOSED SYSTEMS 

Once we have identified a set of statistically interact- 
ing particles, we are not only free to change the energies 
of some or all species to produce new kinds of ordering as 



From Eqs. (21 1 we infer the following expressions for the 



discussed in Sec. IV C we can also impose constraints on 



the number of particles for individual species or combina- 
tion of species. Here we present two simple applications 
that use the particles from Tables [H] and |III| in a closed 
system. Both situations illustrate how effects typically 
attributed to interactions between particles can be at- 
tributed to structures of particles. 



A. Soliton paramagnetism 

We consider a system of N$ solitons that satisfy 
semionic exclusion statistics (Table [TTj) . The solitons pop- 
ulate a lattice of N sites and are able to switch spin ori- 
entation through thermal fluctuations. The solitons have 
spin s± = ±5) therefore, energy e± — ~^\h in a magnetic 
field h. 

Unlike in the ideal Langevin paramagnet, where the 
spins are fixed to lattice sites and thus distinguishable, 
the solitons are indistinguishable floating objects. When 
an individual soliton flips its spin it must shift its position 
left or right by one lattice unit: 



unit 



turn 
mm 



(31) 



In crowded conditions only collective spin-flips of groups 
of solitons are possible. 

We carry out the analysis in the grandcanonical ensem- 
ble and enforce a constant (average) number of particles, 
independent of temperature and magnetic field, via a fu- 
gacity £. Equations (20) in exponentiated form read 



/ -i . m/2 1/2 1/2 



(1 



-)l/2 



H= , (32) 



yielding the physical solutions 



w± = e TH 



sinh 2 H + C~ 2 T sinhiJ 



(33) 



total number of solitons per site and the magnetization 
per site: 



No = 



(N + } + (N-) _ q{(, H) cosh H + sinh 2 H 
~ q((,H)[ q ((,H)+ cosh H}' 



N 



M z = 



(JV+) - (N-) 
2N 



sinh H 



(34) 



(35) 



where q((, H) = \J sinh 2 H + (~ 2 , and the following ex- 
plicit result for the magnetization curves describing soli- 
ton paramagnetism: 



where 

Q(N S ,H) 



M z _ 1 Q(N S , H) + cosh H 
Ns ~ 2Q(N S , H) coth H + sinh H ' 



(36) 



coshi/ 
2N S 



1-N S 



+ J (1 - N s ) 2 + W s tanh^ H . (37) 



Magnetization curves at different soliton densities are 
shown in Fig. [6] 

Corresponding magnetization curves for the s = ^ 
Ising paramagnet, Hamiltonian ([T]) in the limit D — > oo 
previously con side red as an open system of solitons, are 
described by pB) with C" 1 = e 2Kj , Kj = J/4k B T 



and have similar shapes. The common (weak-coupling) 
Langevin limit, M z — | tanh H , is realized for Ns, Kj — > 
and the common strong-coupling limit, M z — ^sgnH, 
for N s -> 1, Kj oo. 

Away from the Langevin limit (dashed curves), the 
magnetization is enhanced due to suppressed thermal 
fluctuations. However, different physical mechanisms are 
responsible in the two cases. In the Ising paramagnet 



0.5 F 




FIG. 6: Magnetization curves of the soliton paramagnet at 
various soliton densities and of the Ising paramagnet at vari- 
ous coupling strengths. 
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individual spin flips are suppressed by a dynamic effect 
(energy barrier). In the soliton paramagnet, individual 
spin flips are suppressed by a kinematic constraint (space 
limitation). 



B. Polymerization 

Consider a system of Nm molecules from a single 
species with a propensity for linking up into polymers. 
The molecules inhabit a one-dimensional lattice of N 
sites. Individual molecules have a certain size and give 
up a fraction of that size with each link to a neighbor. 
Associated with each link is a binding energy e^. Posi- 
tive (negative) values enhance (suppress) polymerization. 
Crowding may drive polymerization even for negative 
binding energies. 

What fraction of molecules are bound in polymers and 
what is the average length of a polymer at any given tem- 
perature? With the host and tag particles from Table III 
we can answer these questions in purely kinematic terms 

Hosts (tit) are monomers or seeds of polymers. Tags 
(44-) are interior molecules of polymers. Each poly- 
mer thus consists of one host and at least one tag 
(tit + 44-=t44-t)- The number of bonds in a polymer 
is equal to the number of tags it contains. The centers 
of two nearest-neighbor monomers are two lattice units 
apart. That distance is reduced to one lattice unit when 
monomers form a dimer |16j . 

The binding energy 6b is accounted for by assigning 
different energies to hosts and tags: €b = £h — £t- Equa- 
tions ( 20 1 for this situation become 



1 = w 2 H {l + w T ) 

c w T (i+w H y 



,-Kb 



w H (l + wt) 
1 + wh 



(38) 



where Kb = eB/ksT. We have absorbed the host energy 
e H in the fugacity (. The physically relevant solutions for 
an open system are 



Wh = WTe KB = — 



1 



We B + \/(We A ' B ) 2 +4C . 

(39) 



The population densities (N m ) = (N m ) /N, m = H,T, 
inferred from Eqs. \21\ are 



(N H ) = w T (N T ) = "T^-, • (40) 

wt{wh + 2) + 1 

Now we switch to a closed system by imposing the con- 
straint (Nh) + (Nt) — Nm- Replacing £ by Nm as 
independent variable we write 



1 



W H 



1 - 2N, 



A I 



2N M 

+ J(l- 2N M ) 2 + 4JVAf (1 - N M )e K * 



(41) 



In Fig. [7] we show the dependence on the molecular 
density Nm of the fraction of molecules that are part of 




FIG. 7: Average length of polymer (main plot) and fraction of 
molecules that are bound in a polymer (inset) versus molec- 
ular population density at various values of en/hsT. The 
dashed lines are for €b = 0. 



a polymer (inset) and the mean polymer length (main 
plot) for various fixed values of scaled coupling. In a 
macroscopic system, N^,N ^> 1, two quantities are 



(N 7 



(N H ) + (N 7 



F P N 



M ■ 



(42) 



For zero binding energy the fraction of polymerized 
molecules increases linearly with Nm and the average 
length of the polymers increases quadratically. Positive 
(negative) binding energies enhance (suppress) this rate 
visibly in both quantities. Large positive binding energies 
lures all available molecules into polymerization and the 
average length increases linearly. Large negative binding 
energy deters all molecules present from polymerization if 
there is enough space for monomers. The average length 
stays near zero, then rises linearly when polymerization 
is forced through crowding. 



VI. CONCLUSIONS 

We have explored several themes of fractional exclusion 
statistics as applied to particles with shapes. The context 
has been Ising spin chains, but the conclusions are readily 
transcribed to different applications. 

Consistency in the combinatorics and statistical me- 
chanics of indistinguishable particles at different levels of 
description has been the main theme. Concerns about 
this issue were raised previously [TUl E] ■ Anghel derived 
a set of necessary conditions that must be satisfied when 
species of statistically interacting particles are merged or 
split [S] . We have identified two distinct types of mergers 
that guarantee consistency in the combinatorial analysis, 
both of which confirm Anghel's rules. Our results confirm 
that Wu's method of statistical mechanical analysis [14] 
provides consistent exact results at each level of descrip- 
tion, i.e. before and after mergers. We have shown evi- 
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dence that any sufficient conditions for splitting species 
must be model specific, in agreement with Refs [9HTT]. 

One subsidiary theme has been about open and closed 
systems of statistically interacting particles. In the Ising 
context the numbers of particles from each species fluc- 
tuate independently. We have described two situations 
where such fluctuations are constrained by a conserva- 
tion law. A second subsidiary theme has been about the 
transcription of motifs for statistically interacting parti- 
cles between models for different physical situations. 

The statistical mechanics of particles with shapes is a 
field wide open to applications of strong current interest 
including granular matter jammed in narrow channels 
[T7] and DNA transformed under tension JT5J, Q15] . 



Appendix A: Ising chain with dimerized coupling 



The s — \ Ising chain with alternating bond strengths, 



2JY 



(Al) 



where Jo = J P + J m is the coupling for bonds with odd I 
and Je = Jp~ Jm for bonds with even I, has four distinct 
T = phases: 



d>_ + 



I tin •••>!, i mt • ■ ■ )i at j > o, j e > o ; 

|mtmt--">2 at J > 0, J E < 0, 
|tmtm--->2 at J < 0, J E > 0, 
| ft U"-)i at Jo <0,Je <0. (A2) 

The periodicity takes the doubling of the unit cell into 
account. One set of particles that generates the spectrum 
of ( Al I on a lattice with 2N sites is closely related with 
the set of particles from Tablc[l]for ([I]) on a lattice with N 
sites. The Hilbert space of both systems has dimension- 
ality 4 N . We transform the inhomogeneous system into a 
homogeneous system by introducing compound site vari- 
ables constructed from pairs of site variables coupled by 
Jo- The four compound site variables combine to 16 
bonds as shown in Table [V] The energy of each com- 
pound bond has two parts: one unit of Jg-bond energy 
(counted fully) and two units of Jo-bond energy (each 
counted half). Each compound site and each compound 
bond is also assigned a spin. 

We now select the physical vacuum at Jo < 0, Je > 0, 

i.e. phase $ |_, as the pseudo- vacuum of a set of particles 

that generate the spectrum of (Al). We take advantage 
of the isomorphism, 



(ft, ft, t, I) 



(ft, ft, n, u), 



(A3) 



that transforms the site variables of ([I]) into the com- 
pound site variables of ( Al ). It maps the pseudo- vacuum 
of the particles listed in Table |T| fo r ([I]) into the the 
pseudo- vacuum selected here for <| All) . 



The particles thus constructed for (Al) have different 
energies e m and different spins s m . The specifications 
Am, a m , g m m' remain the same. The energies of par- 
ticle pairs 15 & 16, 1 & 2, 5 & 6, 3 & 4 remain equal, 
making it possible to carry out the first four mergers as 
in Sec. |III| However, the two motifs of particles 7 and 8 
acquire different energies in the transcription, which ap- 
pears to make it necessary to split them into subspecies. 

Given the indeterminate nature of rules ^ in this 
process, it requires that the combinatorial analysis be 
performed from scratch for the new subspecies including 
their statistical interactions with all other species. In the 
situation at hand, this can be avoided. Consider the four 
motifs fftft, fttft, ftlft, I ftlft of particles 7 and 8. In the 
context of |l]) all four have the same energy at ft. = but 
the first two have spin + 1 and the last two spin — | , pro- 
ducing different Zeeman energies. In the context of ( Al ) 



all four have zero spin but the energy of the first and the 
last is different from the energy of the second and the 
third. Can we switch the second and the fourth motif 
such that particles 7 and 8 each have a definite energy 
again? The answer is affirmative. The switch leaves all 
results of the combinatorial analysis invariant provided 
the merger of species 15 & 16 into 15 is in place. We 
thus end up with 13 species of particles. Their specifica- 
tions are listed in Table IVII 

The statistical mechanical analysis of these particles 
proceeds analogous to that reported in Sec. |IVD| The 
particles have the same motifs (albeit transcribed and 
slightly shuffled) but different energies. The same motifs 



TABLE V: Compound sites associated with Jo-bonds, their 
energies, and their spin content (top) Compound bonds asso- 
ciated with Jfi-bonds of compound sites, their energies, and 
their s pin content (bottom). 
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encode different physical structures and produce different 
kinds of ordering in the limit T —¥ 0. Adding a Zeeman 
term to Q poses no serious challenge. It requires that 
the mergers of 3 & 4, 5 & 6, 1 & 2 be reversed and that 



the merger of 15 & 16 be kept. The statistical mechanical 
analysis then proceeds with 16 species and with Zeeman 
energies —s m h added to e m . 



TABLE VI: Specifications of 13 particles that generate the spectrum of (All from the state (n pv — 2), 
tt-W-tt-U- ■■■)l> = | 1T-0-1T • • - >2 after four mergers and two splits as described in the text: motif, species, energy (relative to 
vacuum), spin, capacity constants, and size constants, statistical interaction coefficients. Segments of £ (compund) vacuum 
bonds Tr-U- ^ 4J-TT have energy E = \t(Jo — Je)- The first species is a compact, the last two are hybrids, and all others hosts. 
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